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Molecular dynamics simulations of nanometer-sized lead clusters have been performed using the 
Lim, Ong and Ercolessi glue potential (Surf. Sci. 269/270, 1109 (1992)). The binding energies of 
clusters forming crystalline (fee), decahedron and icosahedron structures are compared, showing that 
fee cuboctahedra are the most energetically favoured of these polyhedral model structures. However, 
simulations of the freezing of liquid droplets produced a characteristic form of surface-reconstructed 
"shaved" icosahedron, in which atoms are absent at the edges and apexes of the polyhedron. This 
arrangement is energetically favoured for 600-4000 atom clusters. Larger clusters favour crystalline 
structures. Indeed, simulated freezing of a 6525-atom liquid droplet produced an imperfect fee Wulff 
particle, containing a number of parallel stacking faults. The effects of temperature on the preferred 
structure of crystalline clusters below the melting point have been considered. The implications of 
these results for the interpretation of experimental data is discussed. 



I. INTRODUCTION 

It has long been known that atomic clusters of some fee 
metals occur in stable arrangements quite different from 
the bulk crystal structure. Regular non-crystalline struc- 
tures, with £«e-fold axes of symmetry, were first identi- 
fied in_goldulfl, and are known to occur in a variety of 
metalsa. The stability of these arrangements relies on a 
high proportion of surface atoms. For example, a 4 nm 
fee lead cluster has roughly 40 % of its atoms at the sur- 
face. The local environment of surface atoms is quite 
different to those inside, and leads to distinct surface 
properties which play an important role in the overall 
behaviour of the cluster. 

Although the observation of five-fold axes of symme- 
try in gold clusters was unexpected, it was immediately 
recognized that slightly deformed fee tetrahedral units 
could be assembled to form these structurestl. These al- 
ternative structures afford a net gain in binding energy 
for the cluster: the strain of deforming fee tetrahedral 
units, and the energy cost of twinning, is compensated 
for by a higher proportion of atoms in densely packed 
outer fafips than can occur in any faceting of a perfect fee 
particletffl. 

Two characteristic types of non-crystalline structure 
occur: the icosahedron consists of twenty tetrahedra as- 
sembled around a common apex; the decahedron ar- 
ranges five tetrahedra around a common edge. In each, 
deformation of the tetrahedra is necessary to fill the vol- 
ume of the polyhedron. This lowers the symmetry of the 
fee lattice to either rhombohedral (icoj&ahedron) or body- 
centered orthorhombic (decahedron)El. In both cases, 
adjacent tetrahedral faces meet at a twin plane, hence 
the name Multiply Twinned Particle (MTP). Packing 
of atoms to make a given polyhedral form leads to se- 



quences of atom numbers at which .complete outer lay- 
ers, or shells, are "closed" (see e.g£j). Each sequence 
depends on the particular polyhedral model. For exam- 
ple, there are two fee cuboctahedron sequences commonly 
considered: one with equilateral triangular (111) facets, 
the other with hexagonal (111) facets. The closed-shell 
sequences for these are, respectively: 13, 55, 147, 309, 
561, etc; and 38, 201, 586, 1289, etc. The hexagonal 
faceted cuboctahedron is often called the "WulfF' form, 
because, for noble metals, it is close to the so-called Wulff 
polyhedron, which minimizea-|tbc surface energy contri- 
bution in a small fee crystacJL-l In this work we gen- 
erally use the terms Wulff particle and cuboctahedron to 
refer to the hexagonal and triangular forms, respectively. 

Elasticity theory can explain the form and stability of 
MTP structures, however, more detailed theoretical tools 
are needed to explore their properties: molecular dynam- 
ics (MD) has been used for this purpose. Classical MD 
uses a mathematical model for the interaction potential 
between atoms, and the classical equations of motion, to 
simulate the dynamic behaviour of an rt-atom aggregate 
over time under prescribed conditions of energy, temper- 
ature, etc. The technique is computationally intensive, 
but gains in computing power continue to allow more 
elaborate and realistic simulations to be performed. 

Of the fee metals, gold clusters, in particular, have 
been extensively studied by MD using many-body inter- 
atomic potentialsLrtj. Early studies of the energetics of 
relaxed gold clusters suggested that crystalline structures 
were favoured over MTP structurestl. Later investiga- 
tions suggested that a truncated Marks decahedronQ may 
be the preferred low-temperature cluster geometry, when 
there are less than about 250 atoms in-.the cluster, and a 
Wulff-particle form for larger clusterscl. Further studies 
have found that disordered structures prevail for certain 
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small cluster sizesEJ MD has also been used tainvesti- 
gate the coalescence of gold clusters in collisionsEl,|-struc- 
tural transitions in gold clusters prior to meltingc2l and 
the growth of silver clusters from small seed clusters by 
collisior£3. 

The case of lead clusters has received much less at- 
tention than gold. The first comprehensive-jpolecular 
dynamics study, by Lim, Ong and Ercolessio, utilised 
a many-body "glue" potential to compare the energetics 
of certain fee and MTP structures for a range of cluster 
sizes. In particular, these authors compared the ener- 
getics of the cuboctahedron and icosahedron sequences 
of closed-shell structures (the numbers of atoms in the 
icosahedron sequence, conveniently, is the same as those 
in the cuboctahedron sequence). This allowed them to 
directly compare the binding energies at the same num- 
bers of atoms for each structure and they demonstrated 
that the fee crystalline structures (cuboctahedra) were 
favored over icosahedra for all sizes where the glue po- 
tential is expected to hold (considered to be clusters with 
more than 200 atomsLLik However, subsequent simula- 
tions by these authorsta, involving the quenching of a 
large (8217-atom) liquid lead droplet, produced a cluster 
that was "icosahedron-like" . Although known not to be 
the lowest energy arrangement for that number of atoms, 
it was assumed to have formed because of insufficiently 
long equilibration times in the simulation. 

The purpose of this work is to take a more comprehen- 
sive look at lead clusters by considering a wider range 
of structures and more variety in cluster sizes than in 
previous studies. We have also conducted simulations of 
longer duration. We begin by outlining the details of our 
computational procedure. This is followed by a discus- 
sion of the trends in binding energy for crystalline and 
non-crystalline model structures over cluster sizes from 
200 to 6000 atoms. We then consider the melting and 
freezing of fee clusters. Sections [v] and VI then exam- 
ine the structure and energetics of clusters produced by 
freezing liquid droplets. Finally, we look at the character- 
istic diffraction patterns produced by some MD-derived 
cluster structures. 



II. COMPUTATIONAL DETAILS 

A. Molecular dynamics 

Molecular dynamics was performed using a local ver- 
sion of the classical molecular dynamics code ALCMD, 
originally developed by Ames Laboratory. Simulations 
were performed in the canonical ensemble (i.e. constant 
temperature) to produce caloric curves in section IV and 
V. The microcanonical ensemble (i.e. constant energy) 
was used to study re-solidification in section V and VI. 
The time step was chosen as 3.75 fs throughout. Simula- 
tions were typically allowed to equilibrate for 3 x 10 5 time 
steps, or approximately one nanosecond. Longer equili- 



bration times (up to 3 x 10 6 time steps) were used near 
transition points and in re-solidification studies. Fol- 
lowing the equilibration period, quantities reported here 
(such as energies) were obtained by time-averaging over 
1 x 10 5 further steps. 

The inter-atomic potential used is due to Lim, Ong 
and Ercolessii-i This is a many-body glue-type poten- 
tial, given by 
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where 4> is a short-range two-body potential and U is 
a many-body glue term which reflects the effects of the 
conduction electrons in the metal. The rii = ■ p(rij) is 
a generalized coordination where p is some short-ranged 
function. The three functions 4>, p and U have been ob- 
tained by fitting to a number of known properties of lead 
including cohesive energy, surface energy, elastic con- 
stants, phonon frequencies, thermal expansion and melt- 
ing temperatureta. This potential has been used previ- 
ously to model lead clusters£JII3, temperature-depeadent 
surface reconstructions of low-index lead surfacesliZi and 
pre-melting of low-index lead surfacesEEl. While the po- 
tential has been fitted to bulk properties of lead, the clus- 
ters studied here are sufficiently large that (Q) should be 
considered quantitatively reliable. 



B. Common neighbor analysis 

Common neighbor analijsiiEl (CNA) has been used to 
analyse cluster structures&EJ. CNA is a decomposition 
of the radial distribution function (RDF) according to 
the local environment of each pair. We consider that the 
first peak in the RDF represents "bonded" neighbors. As 
such, if r c is the first minimum in the RDF, we classify 
any pair separated by r < r c as a bonded pair. With this 
identification, any pair contributing to the RDF can be 
classified by a set of three indices, ijk, which provide in- 
formation on the local environment of the pair. The first 
index, i, is the number of bonded neighbors common to 
both atoms. The second index, j, is the number of bonds 
between this set of common neighbors. The third index, 
k, is the number of bonds in the longest continuous chain 
formed by the j bonds between common neighbors. Fig- 
ure ^ shows a 421-pair (light grey), with four common 
neighbors (dark grey), two bonds (black) and a longest 
chain of one bond. 
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FIG. 1. Example of a 421-pair (paired atoms in light grey), 
with four common neighbors (dark grey), two bonds (black) 
and a longest chain of one bond. 



Figure [2| shows a 959-atom crystalline structure, 
formed by relaxing a ball of lead atoms "cut" from the fee 
lattice. Noting that the incomplete outer shell has plates 
of atoms on hexagonal (111) faces, we characterise this 
structure as an incomplete Wulff particle. The CNA de- 
composition of the first peak in the RDF for this cluster 
is shown in Figure |3|. Here, the RDF has been averaged 
over 100 time steps after the cluster was brought to ther- 
mal equilibrium at 200 K. The dominant contribution 
to the first peak comes from 421-pairs. This is charac- 
teristic of bulk fee structure: each atom in a fee crystal 
has twelve nearest neighbors and each of the pairs formed 
with these neighbors is a 421-pair. Also apparent are 311- 
pairs, which predominantly arise from atoms on a (111) 
fee face, and 211-pairs, which predominantly arise from 
atoms on a (100) fee face. One other pair type that con- 
tributes only weakly to the first peak (and is thus difficult 
to distinguish in the decomposition shown in Figure ^ is 
the 200-pair, which is associated with atoms at the edges 
of the incomplete outermost atomic layer of this cluster. 




FIG. 2. A crystalline structure formed by relaxing a crys- 
talline ball of 959 lead atoms. The darker atoms denote sur- 
face atoms. Note the incomplete outer shell has plates of 
atoms on hexagonal (111) faces. Hence we characterise this 
structure as an incomplete Wulff particle. 
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FIG. 3. Common neighbor decomposition of first peak in 
the RDF for a 959-atom cuboctahedron (pictured in Figure ^) 
equilibrated at 200 K. 

CNA is useful here because it allows one to distin- 
guish between local atomic arrangements, including fee 
and icosahedral environments. In Table | we have listed 
the classifications of CNA signatures used here to la- 
bel the local environment of an atom (this classification 
is similar but not identical to that used by Cleveland, 
Luedtke and LandmanEll) . We note that these signatures 
are based only on the CNA decomposition of the first 
peak in the RDF. Figure |J shows the results of this clas- 
sification scheme applied to the relaxed crystalline 959- 
atom structure of Figure ^[ Here we see that predomi- 
nantly fec-bulk atoms with (lll)-face atoms, (100)-face 
atoms and (111)/(100) edge atoms are present. 

TABLE I. Description of CNA signatures used. 



? unclassified signature (possibly disordered) 

A fee internal atom 

B fee (111) face atom 

C fee (100) face atom 

D fee (111)/(100) edge atom 

E internal atom at a (111) fee stacking fault 

F internal icosahedral atom (spine or central atom) 

G surface icosahedral apex atom 

H surface icosahedral (111) /(111) edge atom 

I "shaved" surface icosahedral (111) /(111) edge atom 
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FIG. 4. Classification of atoms in the 959-atom incomplete 
Wulff particle (pictured in Figure ^) equilibrated at 200 K. 



III. BINDING ENERGIES OF RELAXED 
GEOMETRIC STRUCTURES 

Lim, Ong and ErcolessiEl compared the energetics 
of the icosahedron sequence with the energetics of the 
cuboctahedron sequence. These structures were relaxed 
to K using the glue potential (Q) and the total energies 
for each pair in this sequence were computed (recall that 
these polyhedra have complete outer shells for 13, 55, 
147, 309, 561, etc atoms). They found that the relaxed 
cuboctahedra were strongly favoured energetically over 
the relaxed icosahedra. 

One cannot directly compare the energetics of struc- 
tures which do not have the same numbers of atoms. In 
general, as cluster size increases, the binding energy per 
atom in the cluster also increases, as the percentage of 
surface atoms decreases. This trend is often described 
by the following relationship between binding energy per 
atom and cluster sizeO: 

E b = A + BN~ 1/3 + CN~ 2/3 . (2) 

By fitting A, B and C to the binding energies of families 
of structures, one can use (|J) to interpolate between clus- 
ters of complete outer shells and thus make comparisons 
between families which do not occur in matching pairs. 

In addition to the cuboctahedron and icosakedron se- 
quences (looked at by Lim, Ong and Ercolessitil) , we also 
consider the sequences of the Ino decahedrao and Wulff 
cuboctahedra, making use of ^ to compare energies by 
interpolation. The Wulff construction minimizes the sur- 
face energy in a small crystal by taking into account the 
relative surface energy of distinct crystal facetscJeJ. For 
this reason, one expects Wulff clusters to be energetically 
favoured with respect to the alternative cuboctahedron 
sequence. However, our results will show that any differ- 
ences in binding energy are too small to be significant. 
We also l&pk at clusters extracted as a sphere from the 
fee lattices. The number of atoms in these ball clusters 



typically does not coincide with that of a closed-shell 
structure. Nonetheless, when relaxed from the bulk they 
tend to resemble the Wulff form (see Figure ||) as noted 
in the previous section; excess atoms in the incomplete 
outer shell contribute to hexagonal (111) facets. 



The energies of the three types of crystalline model 
examined, are shown in Figure ^together with the inter- 
polations based on equation (|^). For each, we have con- 
sidered a sequence of sizes, relaxing atomic configurations 
to K and computing the binding energy per atom. As 
might be expected, the incomplete Wulff particles were 
less favoured than the cuboctahedra or Wulff particles, 
because of the incomplete (111) facets. However there 
appears to be little distinction energetically between the 
Wulff and cuboctahedron structures. This contrasts with 
similar comparisons in Lennard- Jones clusterso and in 
nickel clusterscJ, both of which identify the Wulff parti- 
cles as favorites. 




FIG. 5. Binding energies of relaxed crystalline clusters (□ 
Wulff particles, • cuboctahedra and o incomplete shell Wulff 
particles) relative to the fit to the energy of the cuboctahe- 
dron sequence. For reference, the numbers of atoms of each 
member in the cuboctahedron sequence are shown. 



In Figure |6J, we compare the energy of the MTP struc- 
tures (decahedra and icosahedra) with cuboctahedra and 
incomplete Wulff particles. 
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FIG. 6. Binding energies of relaxed clusters (• cuboctahe- 
dra, A Ino decahedra, o incomplete shell cuboctahedra, O 
icosahedra ) relative to the fit to the energy of the cubocta- 
hedron sequence. For reference, the numbers of atoms of each 
member in sequence of cuboctahedron clusters are shown. 

From Figure |(| we see that cuboctahedra are favoured, 
in general, over MTP structures and fee structures 
with an incomplete outer shell. Ino decahedra are 
also favoured over the incomplete cuboctahedra, and all 
these structures are favoured over icosahedra. In gen- 
eral, these, results support the findings of Lim, Ong and 
ErcolessiliJ, that the cuboctahedra are the most energet- 
ically favoured structures. Nonetheless, for certain clus- 
ter sizes the relaxed Ino decahedra are the lowest en- 
ergy structures considered. This is unexpected, because 
the Ino decahedron is usually regarded as an unfavorable 
atomic arrangement; it is the more complex re-entrant 
faceted model of Marks that is-xonsidered energetically 
competitive in small particlesE3d. However, in the next 
section we see that none of these geometric models con- 
stitute the lowest energy structures for a range of cluster- 
sizes. 



teresis" in the melting/freezing point is typical in MD 
simulation, an effect of similar magnitude has been ob- 
served experimentally in lead clusters, where liquid lead 
clusters have been undp.-cooled by as much as 120 K 
without re-solidificationEa. 
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FIG. 7. Caloric curves for fee ball clusters: 555-atom clus- 
ter, □ ; 959-atom cluster, o ; 1721-atom cluster, A . The 
empty symbols indicate heating and the filled symbols repre- 
sent cooling. 

As clusters were heated towards the melting point, a 
number of interesting structural changes could be iden- 
tified from the CNA classification. Figure || compares 
the CNA results of a 959-atom incomplete Wulff particle 
at temperatures 250-350 K, prior to melting. One can 
identify a general trend towards disorder as the number 
of unclassified atoms (signature ?) increases. However, 
at 350 K we see a sharp increase in atoms with signa- 
ture E, predominantly in the outermost complete shell 
of the crystal. The signature E atoms arise as the in- 
complete (111) facets above these atoms shift to form 
stacking faults at the surface. 



IV. MELTING AND FREEZING OF CLUSTERS 

In a later study, Lim, Ong and ErcolessiEl considered 
the simulated quenching of an 8217-atom liquid clus- 
ter and observed the formation of an imperfect icosa- 
hedron upon solidification at constant energy. Here, we 
have performed a more general investigation of the melt- 
ing and freezing of lead clusters. Figure shows the 
caloric curves for several sizes of incomplete Wulff clus- 
ters (which were first relaxed to their K structures then 
heated). The caloric curves were constructed using a 
sequence of constant temperature simulations. The en- 
ergy varies linearly with temperature on either side of the 
melting/freezing points, where a sharp change in energy 
occurs (latent heat of melting). The melting points of 
the clusters show a general increase with cluster size as 
expected, with all melting points lying below that of the 
bulk melting point (approximately 600 K). While "hys- 
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FIG. 8. Classification based on CNA signatures of atoms 
in the 959-atom cluster for temperatures 250-350 K. 

Returning again to the caloric curves (Figure 0) , we see 
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that, once re-solidified, clusters have a lower energy than 
the original relaxed structures, indicating that atomic re- 
arrangement has probably occurred. For cluster sizes of 
between 555 and 3500 atoms, all the re-solidified struc- 
tures examined were found to be icosahedron-likc. Fig- 
ure shows one example of a re-solidified 959-atom clus- 
ter. An axis of five- fold symmetry is clearly visible at 
the surface and predominantly (111) faces are exposed. 
Figure ^ compares the binding energy of a number of 
re-solidified structures with cuboctahedra. As was ap- 
parent in Figure ^, the re-solidified clusters appear to be 
favoured energetically over cuboctahedra for a range of 
sizes from approximately 600 atoms up to about 4000 
atoms. The interpolations of the binding energy of each 
family, suggest that the fee structures may become more 
favorable above the 4000 atom size. One would expect to 
see a crossing of the two curves at large cluster sizes, sim- 
ply because the bulk structure should be favored in this 
limit. The largest re-solidified structure in this study was 
the 6525-atom cluster shown in Figure 11 (this was pro- 
duced using the constant energy method detailed in sec- 
Here, we clearly see an imperfect fee Wulff par- 
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ticle structure emerging, which is consistent with these 
considerations (we consider the structure further in sec- 
tion Ivjl). 
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FIG. 10. Re-solidified 959-atom cluster. Note the axis of 
five- fold symmetry with "shaved" edges (darker atoms). Also 
visible are exposed excess atom layers on (111) faces. 

Lim, Ong and ErcolessiS report a similar icosahedron- 
like structure that formed upon re-solidification of an 
8217-atom cluster (using the constant energy method 
similar to the method detailed in the next section) . From 
considerations of the binding energies of geometric icosa- 
hedra and cuboctahedra (see Figure ^), they expected 
the quenched liquid droplet to form a fee cuboctahe- 
dron. The appearance of an imperfect icosahedron was 
attributed to an insufficient relaxation time. While the 
interpolation (||) indicates that for clusters of sizes of 
up to 4000 atoms, the icosahedron-like structures are 
favoured over fee, extrapolating (||) suggests that clus- 
ters above 4000 atoms favor crystalline structures. We 
conclude that the failure of earlier studies to find fee 
structures on re-solidification is in part due to shorter 
relaxation times. We note that the relaxation times used 
here, are-several times longer than those of Lim, Ong and 
Ercolessilij. 
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FIG. 9. Comparison of binding energies for cuboctahedra 
(•) and re-solidified icosahedron-like structures (v/) relative 
to the fit to the energy of the cuboctahedron sequence. The 
relative binding energies of icosahedron-like clusters for sizes 
555, 923, 959, 1721, 2057 and 3428 are shown here. Apart 
from the 555-atom cluster, these structures are energetically 
favored over cuboctahedra. 




FIG. 11. Re-solidified 6525-atom cluster showing faceting 
typical of the Wulff particle, a cuboctahedron with hexagonal 
(111) faces. 

The re-solidified icosahedra in this section were pro- 
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duced using constant temperature simulations. In the 
following section, we show that these icosahedra are also 
produced from liquid droplets by freezing at constant en- 
ergy. We will also characterise the structure of these 
icosahedra using CNA analysis. 



V. RE-SOLIDIFIED ICOSAHEDRA 



A. Structure of re-solidified icosahedra 



In this section, we examine the structure of a "shaved" 
re-solidified icosahedron cluster produced using constant- 
energy re-solidification. Starting with a perturbed 2057- 
atom Wulff particle (each atom was perturbed in a ran- 
dom direction with mean displacement 0.1 A), the clus- 
ter was melted and then brought to equilibrium at 900 K 
for approximately 1 ns. The cluster was then quenched 
rapidly below the melting point, where it was allowed to 
equilibrate at constant energy for approximately 10 ns. 
The structure was then slowly relaxed to K where the 
binding energy was computed and the structure was de- 
termined with CNA analysis. The re-solidification proce.- 
dure is similar to that used by Lim, Ong and ErcolessitS. 

Figure |l2| compares the CNA classification of a 2057- 
atom icosahedron, before melting, with a re-solidified 
2057-atom cluster. We see that disorder seems to have in- 
creased with re-solidification (signature ?), and the num- 
ber of internal fee atoms (signature A) and (111) face 
atoms (signature B) have also decreased. Also of note 
is the appearance of (lOO)-faceting (signature C) as well 
as the removal of edge atoms along a twin plane (lead- 
ing to signature I atoms). An increase in the number of 
hep- type (signature E) atoms is also apparent. 

This re-solidified icosahedron is a surface- 
reconstructed, or "shaved", icosahedron, where the 
(111)/ (111) facet edges and the icosahedral apex atoms 
on the outer shell have been redistributed in the cluster, 
and the (111) layers have moved to form stacking faults 
with respect to the internal fee tetrahedra of the icosa- 
hedron (leading to an increase in signature E atoms). 
While the disappearance of the interior fee atoms ap- 
pears, in part, to be due an increase in disorder, the 
shaving of (111) / (111) edges gives the internal fee atoms 
below these edges icosahedral coordination, resulting in 
these atoms being classified as "extra" internal icosa- 
hedral atoms (signature F). This "shaved" icosahedron 
surface is the same surface reconstruction as that noted 
by Doye and Wales in theic, investigation of icosahe- 
dra using a Morse potentialEa. Their study found that 
geometric icosahedra underwent thermally-induced sur- 
face reconstructions prior to melting to form the same 
"shaved" icosahedra seen here. 
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FIG. 12. The CNA classification of a 2057-atom icosa- 
hedron compared with the classification of a 2057-atom 
re-solidified cluster. 

The average binding energies of the atoms according to 
classification is given in Table 0. The replacement of the 
(111)/ (111) facet edges, with an average binding energy 
of 1.463 eV/atom, by shaved facet edges with an average 
binding energy of 1.631 eV/atom, considerably increases 
the total binding energy. This compensates for the in- 
crease in internal icosahedral atoms (signature F, 1.997 
eV/atom) at the expense of internal fee atoms (signa- 
ture A, 2.016 eV/atom). This rearrangement of surface 
atoms, with an icosahedral internal structure, is ener- 
getically favoured over the crystalline cuboctahedra for 
cluster sizes from 600 to 4000 atoms (Figure ^) . In their 
original study of lead clusters using the potential (Q), 
Lim, Ong and ErcolessiEJ concluded that fee crystalline 
structures were energetically favoured over this range. 
However, they did not consider the effect of the "shaved" 
surface reconstruction which, as is evident here, stabilises 
the icosahedra. 

TABLE II. Comparison of average binding energies 
of atoms in an icosahedron prior to melting and after 
re-solidification. 





Cluster 






Closed-shell 


Re-solidified 


Signature 


N 


E 


N 


E 






(eV/atom) 




(eV/atom) 


? surface 






169 


1.541 


? interior 






142 


1.975 


A 


700 


2.016 


437 


2.007 


B 


420 


1.625 


165 


1.703 


C 






72 


1.638 


D 






8 


1.471 


E 


630 


2.010 


645 


2.010 


F 


85 


1.992 


188 


1.997 


G 


12 


1.056 


10 


1.674 


H 


210 


1.463 


7 


1.533 


I 






214 


1.631 
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B. Re-solidification of liquid droplets 



C. Structural instability of geometric icosahedra 



In this section we have examined a total of twenty- 
five 2057-atom cluster structures produced by re- 
solidification at constant-energy using the scheme out- 
lined in the previous section. This enables us to in- 
vestigate the reproducibility of the "shaved" icosahe- 
dron structures. The distribution of binding energies 
for the twenty-five re-solidified clusters is shown in Fig- 
ure [l3]. Two re-solidified clusters were characterised 
by CNA analysis as Wulff particles, with energies of 
1.882 eV/atom (equal to the binding energy for a per- 
fect geometric Wulff particle). The remaining twenty- 
three were "shaved" icosahedra (energies above 1.884 
eV/atom). The two "shaved" icosahedra with energy 
1.884-1.885 eV/atom had surface defects (one edge not 
shaved). However, the majority of the "shaved" icosa- 
hedra are clustered about an average energy of 1.887 
eV/atom. Thus we can see that with equilibration times 
of 10 ns kinetic trapping in less optimal structures is un- 
common. 
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FIG. 13. Distribution of binding energies for 25 re-solidified 
2057 atom clusters. Two re-solidified clusters were identified 
as Wulff particles (energies of 1.882 eV/atom) while the re- 
maining 23 were identified as "shaved" icosahedra (energies 
above 1.884 eV/atom). 

The repeated occurrence of the "shaved" icosahedron 
structure in these trials suggests that it is the preferred 
structure at this size for the potential (0). This method 
did not produce any structures which were energetically 
more favoured than the "shaved" icosahedra in any of the 
trials. Furthermore, these "shaved" icosahedra appear 
upon re-solidification using both constant temperature 
(section IV) and constant energy (section V) simulations 
at a variety of temperatures, cooling rates and equilibra- 
tion times (up to the 10 ns used in this section). In the 
next section, we see that by heating geometric icosahe- 
dra, a transformation to "shaved" icosahedra occurs prior 
to melting. 



In this section we show that geometric icosahedra can 
be induced by heating (below the melting point) to form 
"shaved" icosahedra in the same manner as that seen by 
Doye and Waleso. Figure [li] shows the caloric curve for 
two icosahedron clusters of 923 and 2057 atoms. These 
clusters have been prepared by first relaxing a geometric 
icosahedron to K. The clusters were then brought to 
equilibrium at 100 K, and then heated in steps of 25 K 
to 300 K. At each temperature, the clusters were allowed 
to equilibrate over 5 x 10 5 time steps at constant temper- 
ature. Between 200 K and 250 K, the energy of clusters 
is seen to drop sharply, suggesting a change in cluster 
structure. 




FIG. 14. Caloric curve for 923 (□) and 2057 (•) atom icosa- 
hedron clusters. The 923-atom cluster undergoes a structural 
transition at 250 K. The 2057-atom cluster undergoes a struc- 
tural transition at approximately 225 K. 

The new structure that appears at temperatures above 
this transition point is essentially the same as the 
"shaved" icosahedron clusters obtained in the previous 
section after re-solidification. We see that geometric 
icosahedra relax to the equilibrium configuration (the 
"shaved" icosahedron) once heated to a certain temper- 
ature. The temperature at which this relaxation occurs 
is seen to depend on cluster size (which determines the 
energy barrier for this process) but will also depend on 
the time scale of the simulation. 



VI. RE-SOLIDIFIED WULFF PARTICLE 

Here we discuss the structure of the re-solidified 6525- 
atom crystalline structure (shown in Figure pr| ) produced 
by freezing a 6525-atom icosahed ron u sing the constant 
energy method outlined in section V A . Apparent in Fig- 
ure [y] is the appearance of hexagonal (111) faces. This 
is interesting because a 6525-atom cluster can form a 
perfect, closed-shell, cuboctahedron with triangular faces 
but not a perfect Wulff particle (cuboctahedron with 
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hexagonal faces). The average binding energy per atom 
of this re-solidified cluster is 1.92953 eV, as opposed to 
the binding energy of 1.92982 eV for a model 6525-atom 
cuboctahedron with triangular (111) faces. The differ- 
ence of 0.3 meV per atom, in favour of the cuboctahedron 
with triangular (111) faces, does not explain the prefer- 
ence for the Wulff form. However, it is probable that 
the Wulff particle represents a local minimum in the en- 
ergy in configuration space and that the liquid droplet is 
closer in configuration space, due to spherical symmetry, 
to this form of cuboctahedron. A lengthier relaxation 
period upon re-solidification may see the cluster form a 
cuboctahedron with triangular (111) faces. 




The CNA classification of the re-solidified cluster is 
given in Figure The profile is essentially fee, ex- 
cept for the large number of type-E atoms. These arise 
from internal stacking-fault defects in the structure, ei- 
ther when a twin plane occurs (the fee (111) pack- 
ing sequence goes, e.g.: ...ABCACBAC...), or when a 
skip (deformation fault) in the sequence occurs (e.g.: 
...ABCBCABC...). Note that the stacking fault planes 
are all oriented in the same direction. The signature 
E classification (Table is based on the decomposition 
of the first peak of the RDF, and does not distinguish 
between the two. types of defect (although it would be 
useful to do soE!l). The E-type fault planes are shown 
in Figure [H3|. Ignoring the cap regions, there is a set 
of four consecutive fault planes in the upper half of the 
cluster which is a series of deformation faults that has 
effectively produced a narrow band of hexagonal-close- 
packed structure. In the lower half of the cluster, a pair 
of deformation fault planes occur, and below them is a 
single twin-plane fault. 
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FIG. 16. Atoms which show signature E (stacking fault or 
twinning plane) in the re-solidified 6525-atom cluster. 



VII. DIFFRACTION PATTERNS 



Diffraction is a valuable experimental probe of 
nanoparticle structureEj. The diffraction pattern of a 
nanoparticle sample is a "powder" pattern with radially 
symmetric intensity that can readily be calculated from 
the radial distribution function, using the Debye equa- 
tion of kinematic diffraction theoryE3. The so-called "in- 
terference function" is closely related to the distribution 
of intensity in an experimental diffraction pattern and is 
given by 



J(s) 



N ^ 



sin(27rsr m „) 

2 7T S7*777,7i 



(3) 



FIG. 15. The CNA classification of a 6525-atom 
re-solidified cluster. 



where the scattering parameter s = 2sin(0)/A, with 9 
equal to half the scattering angle and A the wavelength 
of the radiation. N is the number of atoms in the cluster 
and r mn is the inter-atomic distance between atoms m 
and n. J(s) is dimcnsionless; it is related to the distribu- 
tion scattered intensity, but is independent of the type of 
radiation and- for monoatomic clusters, the atomic scat- 
tering factorES. 

It is generally difficult to interpret diffraction measure- 
ments on nanoparticles because of the likely occurrence of 
non-crystalline and imperfect structures, as well as sam- 
ple size distribution and inhomogeneity, so the standard 
techniques of x-ray crystallography do not apply. Deter- 
mination of structure is usually based on the careful com- 
parison of experimental data with diffraction patterns 
calculated from models of particle structures. When a 
close correspondence between observations and model- 
based profiles is obtained, the principal characteristics 
of the model structure can generally be inferred to ex- 
ist in the actual cluster sample. The best known and 
most successful example of this approach is probably the 
electron diffraction studies ofr*ase-gas clusters that were 
interpreted using MD modelsE2lE-X 
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Typically in experimental studies, geometric models 
are used to provide the needed candidate diffraction pro- 
files. Thus it is of interest to consider how MD mod- 
els of structure compare with geometrical models. With 
this in mind, we consider the interference functions of 
the two larger clusters that have been generated by re- 
solidification in this study. Figure [jj] shows the pat- 
tern obtained from the 3428-atom cluster, and Figure |l8| 
shows the pattern for the 6525-atom cluster. 

The 3428-atom cluster was obtained after melting a 
decahedron, so in Figure [l?] we also show the interfer- 
ence pattern of a 9-shell decahedron (3428 atoms - the 
initial structure) and a 9-shell icosahedron (2869 atoms). 
Several comments are in order. The MD simulated clus- 
ter has a smoother interference function that either of 
the geometric models; the features are broader and the 
MD profile is also noticeably attenuated with increasing 
s. This is all consistent with a high degree of static dis- 
order in the MD structure. The four significant peaked 
features in the MD curve (at: s ~ 0.37, 0.6, 0.7, 0.9) 
are more closely matched by the profile of the icosahe- 
dron, although the decahedron also has peaks at these 
positions. In either case, the position of the features in 
the MD profile is shifted to higher s (the shift is about 

0.01 A ), indicating a contraction of the length scale for 
the MD structure of the order of 3 %, compared to the 
geometric forms. 

Clearly, from the viewpoint of a diffraction measure- 
ments, our MD cluster would be better matched by 
icosahedral models than either fee structures or decahe- 
dra. Diffraction therefore supports CNA in this case, 
although they are fundamentally different types of anal- 
ysis: CNA characterises the local environment of atom 
pairs, whereas diffraction is a reciprocal space view of the 
radial distribution function, and hence tends to identify 
longer range correlations in a structure. 
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FIG. 17. Interference function for the 3428 atom 
re-solidified cluster ( — ) compared to that for the original dec- 
ahedron (- ■ -) and a 2869 atom icosahedron ( ). 

The 6525-atom cluster has been identified by CNA as 
having the fee structure of the bulk. Nevertheless, Fig- 
ure hq shows, again, the attenuation of the structure's 
interference function with increasing s, indicating sub- 
stantial disorder. Figure |l8| also shows an interference 
function calculated for an ensemble of perfect twinned 
structures. These were fee ball clusters with random 
stacking faults along the (111) axis. A set of 25 interfer- 
ence functions from such models were averaged to pro- 
duce the profile shown. The main features of the MD 
structure interference function coincide with the model 
ensemble pattern, although the relative intensities differ 
substantially. It is also interesting that the location of 
features in the MD profile and the static models profile 
are not shifted, as they were for the 3428-atom struc- 
ture. This means that the lattice parameter predicted 
by the MD simulation agrees with the bulk value used to 
construct the models. 




0.50 1.00 . 1.50 

Scattering parameter (A 1 ) 

FIG. 18. Interference function for the 6525 atom 
re-solidified cluster ( — ) compared to that of an ensemble of 
twinned structures ( ). 



VIII. DISCUSSION 

We have considered the binding energies of relaxed lead 
clusters for a range of particle sizes and cluster struc- 
tures using the "glue" potential ([!]). For medium clus- 
ter sizes (600 to 4000 atoms), we found that a surface- 
reconstructed icosahedron structure was the lowest en- 
ergy structure of those considered. These are shaved 
icosahedra, with atoms removed from the edges and miss- 
ing five- fold apex atoms, and surface atoms relocated 
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to form a stacking fault with respect to the internal foe 
tetrahedra of the icosahedron. Earlier theoretical studies 
of the energetics of MTPs by MarksJ, suggested that the 
removal of icosahedral apex atoms would lead to improve- 
ments in energy. This was supported by_atomistic energy 
calculations far Lennard-Jones clusterscJ, but not in the 
case of nickelc3. In practice, it would be virtually impos- 
sible to observe this shaving by microscopy, as explained 
by Buffat, Flueli, Spycher, Stadelman and BorelE3, and 
we are not aware of any reported observations. However, 
when indeed microscopy can identify cluster structure.-as 
icosahedral, particles generally appear to be imperfecta 
This raises the interesting possibility that the MD struc- 
tures grown here may compare well with typical high- 
resolution images of some small particles. 

In this work we did not identify signs of icosahedral 
precursors to melting as reported in simulations of gold 
clustersL3. These, however, might be revealed by con- 
stant energy simulations (as used in the studies of goldEj) 
of fee clusters near the melting point. We speculate here 
that these precursors, if found, will tend to resemble our 
shaved icosahedral structure more clo sely than the geo- 
metric icosahedra. Indeed, in section VC, we saw that 
by heating, geometric icosahedra could be induced to un- 
dergo a surface reconstruction to the shaved icosahedra 
prior to melting. 

Our finding that icosahedron, rather than fee, struc- 
tures are preferred for certain sizes contrasts with an ear- 
lier MD study of lead that used exactly the &am£ poten- 



tial model, but with shorter simulation timesEj,E3. While 



at the time of that study there were no experimental re- 
ports of lead icosahedra, a recent investigation has indi- 
cated that they do occur abundantlyLj. Electron diffrac- 
tion experiments on a beam unsupported lead clusters 
has identified both decahedral and icosahedral signatures 
in the diffraction patterns for beams with average cluster 
size ranging between approximately 3 and 6 nm. Pre- 
liminary analysis of those results has relied on geomet- 
ric models of cluster structures, however the icosahedral 
component in the data is strong and unambiguous. On 
the other hand, the decahedral component, identified in 
larger clusters, may also be compatible with the kind of 
imperfect and faulted structure that has been obtained 
in the present study for the 6525-atom cluster. 

The imperfect structure of this 6525-atom cluster is 
interesting from the perspective of experimental obser- 
vations. The cluster's main feature is a series of (111) 
stacking faults, all parallel to the same axis. These are 
mainly deformation faults, although one twin plane is 
present. Very recently, de Feraudy and Torchet com- 
mented that parallel stacking faults were most unlikely 
to occur in unsupported argon clustersEj. Our findings 
suggest that, for lead, this may not be the case. If so, this 
poses yet more difficulty in interpreting cluster diffrac- 
tion measurements. De Feraudy and Torchet point out 
that if faulted clusters of this type are to be consid- 
ered, then crystallographic-based calculations based on 
a faulted bulk crystal model are not satisfactory: de- 



tailed atomistic modeling is required0. However, our 
simple attempt to produce approximate atomistic mod- 
els of clusters containing parallel fault planes have been 
unsatisfactory. As shown in Figure [l^, the degree of ac- 
cord with the MD structure is quite poor, and would not 
be accepted in an objective analysis of real experimental 
data. 

In conclusion, MD simulations of lead clusters have 
identified new non-crystalline low energy structures after 
re-solidification of medium-sized liquid droplets. These 
forms are surface reconstructed icosahedra; their bind- 
ing energies are higher than the corresponding perfect 
(relaxed) icosahedron as well as the Wulff form. Wc 
have also shown that sufficiently large liquid droplets re- 
solidify into imperfect crystalline structures, resembling 
the Wulff particle, provided adequate equilibration times 
are used. We have made a preliminary investigation of 
the effects of temperature on structure, noting the ap- 
pearance of stacking faults at the surface of incomplete 
cuboctahedra and a similar solid-to-solid structural tran- 
sition in geometric icosahedra. 



ACKNOWLEDGMENTS 

We would like to acknowledge the use of the Ames 
Laboratory Classical MD code. We are also grateful to 
F. Ercolessi for making the lead "glue" potential freely 
available. We would also like to thank Jonathon Doye 
for his comments. 



1 S. Ino, J. Phys. Soc. Jpn 21, 346 (1966). 

2 J. G. Allpress and J. V. Sanders, Surf. Sci. 7, 1 (1967). 

3 L. D. Marks, Rep. Prog. Phys. 57, 603 (1994). 

4 S. Ino, J. Phys. Soc. Jpn 27, 941 (1969). 

5 L. D. Marks, Phil. Mag. A 49, 81 (1984). 

6 C. Y. Yang, J. Cryst. Growth 47, 274 (1979). 

7 F. Ercolessi, W. Andreoni, and E. Tosatti, Physical Review 
Letters 66, 911 (1991). 

8 L. J. Lewis, P. Jensen, and J.-L. Barrat, Physical Review 
B 56, 2248 (1997). 

9 C. L. Cleveland, U. Landman, T. Schaaff, and M. N. 
Shafigullin, Physical Review Letters 79, 1873 (1997). 

10 C. L. Cleveland, W. D. Luedtke, and U. Landman, Physical 
Review Letters 81, 2036 (1998). 

11 K. Michaelian, N. Rendon and I. L. Garzon, Physical Re- 
view B 60, 2000 (1999). 

12 J. M. Soler, M. R. Beltran, K. Michaelian, I. L. Garzon, 
P. Ordejon, D. Sanchez-Portal and E. Artacho, Physical 
Review B 61, 5771 (2000). 

13 F. Baletto, C. Mottet, and R. Ferrando, Physical Review 
Letters 84, 5544 (2000). 

14 H. S. Lim, C. K. Ong, and F. Ercolessi, Surface Science 
269/270, 1109 (1992). 



11 



15 T. P. Martin, Phys. Rep. 273, 199 (1996). 

16 H. S. Lim, C. K. Ong, and F. Ercolessi, Computational 
Materials Science 2, 495 (1994). 

17 C. P. Toh, C. K. Ong, and F. Ercolessi, Physical Review B 
50, 17507 (1994). 

18 A. Landa, P. Wynblatt, H. Hakkinen, R. N. Barnett and 
U. Landman Physical Review B 51, 10972 (1995). 

19 F. Ercolessi, M. Parrinello, and E. Tosatti, Philosophical 
Magazine A 58, 213 (1988). 

20 A. S. Clarke and H. Jonsson, Physical Review E 47, 3975 
(1994). 

21 C. L. Cleveland, W. D. Luedtke, and U. Landman, Physical 
Review B 60, 5065 (1999). 

22 C. Herring, Phys. Rev. 82, 87 (1951). 

23 C. Herring, in Structure and Properties of Solid Surfaces, 
edited by R. Comer and C. S. Smith (The University of 
Chicago Press, Chicago, 1955), Chap. 1. 

24 R. Defay and A. Prigigine, Surface Tension and Adsorption 
(Longmans, Green and Co., London, 1966), Chap. 17. 

25 B. Raoult, J. Farges, M. F. De Feraudy, and G. Torchet, 
Phil. Mag. B 60, 881 (1989). 

26 C. L. Cleveland and U. Landman, J. Chem. Phys. 94, 7376 



(1991). 

27 L. D. Marks, J. Cryst. Growth 61, 556 (1983). 

28 K. F. Peters, J. B. Cohen, and Y.-W. Chung, Phys. Rev. 
B 57, 13430 (1998). 

29 J. P. K Doye and D. J. Wales, Zeit. Phys. D 40, 466, (1997). 

30 B. E. Warren, X-ray diffraction (Dover, New York, 1990). 

31 D. Zanchet, B. D. Hall, and D. Ugarte, Characterisation of 
nanophase materials (Wiley- VCH, Weinheim, 2000), pp. 
13-36. 

32 A. Guinier, X-ray diffraction in crystals, imperfect crystals 
and amorphous bodies (Dover, New York, 1994). 

33 J. Farges, M. F. de Feraudy, B. Raoult, and G. Torchet, J. 
Chem. Phys. 78, 5067 (1983). 

34 J. Farges, M. F. de Feraudy, B. Raoult, and G. Torchet, J. 
Chem. Phys. 83, 3491 (1986). 

35 P. Buffat, M. Fliieli, R. Spycher, P. Stadelman and J-P. 
Borel, Faraday Discussions of the Chemical Society 92, 173 
(1991). 

36 M. Hyslop, A. Wurl, S. A. Brown, B. D. Hall, R. Monot, 
Eur. Phys. J., in press. 

37 M.-F. de Feraudy and G. Torchet, J. Cryst. Growth 217, 
449 (2000). 



12 



